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Abstract 

We give a geometrical derivation of the Dirac equation by considering a spin-^ particle travelling 
with the speed of light in a cubic spacetime lattice. The mass of the particle acts to flip the 
multi-component wavefunction at the lattice sites. Starting with a difference equation for the 
case of one spatial and one time dimensions, we generalize the approach to higher dimensions. 
Interactions with external electromagnetic and gravitational fields are also considered. One logical 
interpretation is that only at the lattice sites is the spin A particle aware of its mass and the 
presence of external fields. 
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I. INTRODUCTION 


There are different ways to derive the Dirac equation. But probably there is no derivation 


more elegant than the one Dirac gave in his book.[l] The derivation based on Wigner’s 
analysis of the irreducible unitary representation of the “Poincare group” (the covering group 
of the inhomogeneous proper othochronous Lorentz group) certainly is also important. 2] 
There is another intriguing derivation which Feynman gave (for the SU^l, 1) case with one 
spatial dimension and pne time dimension) in his class 1 and which was given as a problem 
in his book with Hibbsjsj]. 

In this paper we give another derivation of the Dirac equation. Our approach bears some 
resemblance to Feynman’s and is based on Dirac’s observation that the instantaneous ve¬ 
locity operators of the spin-| particle (hereafter called by the generic name “the electron”) 
have eigenvalues ±c and that they anticommute. 2 (Hereafter, unless clarity demands oth¬ 
erwise, we set the speed of light c, as well as Planck’s constant h, equal to unity.) We 
assume spacetime to be “filled” with a four-dimensional cubic lattice with lattice length 
Aa; = Ay = A z — At — l. While it is natural to take l to be the Planck length (~ 10~ 33 
cm), we will simply take it to be a length very small compared to the electron’s Compton 
wavelength, the only intrinsic length available in the problem for a free electron. For the 
resulting difference equation, the zeroth order term in At gives a trivial identity while the 
first order term yields the Dirac equation. 

It is interesting that the Dirac equation is invariant under rotations and Lorentz trans¬ 
formations, while the underlying spacetime lattice is not. This situation is one which is not 
unfamiliar in the spatial dimensions in condensed matter physics. (But a related and more 


intriguing result was that found by Snyder more than half a century ago, who showed that 
spacetime being a continuum is not required by Lorentz invariance.) In our approach, an 
electron’s propagation through spacetime can be visualized as consisting of two steps: the 


1 L. Brown, private communication (2001). 

2 The Hamiltonian for a free electron is given by H = caipt + p 3 mc 2 (in Dirac’s notation) with anticom¬ 
muting on and p 3 . The i-th component of the velocity is ^ H] = cat- Thus ^ has as 

eigenvalues ±c , corresponding to the eigenvalues ±1 of cm This result is actually implied by the uncer¬ 
tainty principle. Dirac[lj] also shows that cm consists of two parts, a constant part c 2 piH -1 , connected 
with the momentum by the classical relativistic formula, and an oscillatory part whose frequency is high, 
being at least me 2 /ttTi. 
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electron transfers from one spatial lattice site to a neighboring site in one unit of time (thus 
travelling with the speed of light) and at a lattice site the electron (its multi-component 
wavefunction, to be more precise) is “flipped” by a mass operator 3 and interacts with exter¬ 
nal fields. 

In the next section, we consider the case of one spatial dimension (and one time dimen¬ 
sion). We treat the case of higher spatial dimensions in Section III. Interactions with external 
electromagnetic and gravitational fields are considered in Section IV. Further discussions are 
given in the last Section. 


II. 50(1,1): (1 + l)-DIMENSIONAL CASE 


We assume that the electron of mass m moves with the speed of light from one lattice site 
to a neighboring (spatially left or right) site with time t always increasing on the “cubic” 
spacetime (z,t) lattice. The wavefunction has two components 


z,t) 




( 1 ) 


where denotes the component arriving from the event (z — At, t — At) while means 
arriving from (z + At, t — At). 

Next we assume that, at the lattice site ( z,t ), the arriving components are partially 
turned around by a unitary matrix: 


' Mz,t) \ 

/ 

= F 

K *p-(z,t) J 

V 


ip+(z — At, t — At) 
^-(z + At, t — At) 


( 2 ) 


with the “flip operator” F defined by 


F = e 


—iFm.At. 


( 3 ) 


Here F is a hermitian 2x2 matrix which we call the “flip matrix” and give the most obvious 
form 


F — (J x = (J i — 



(4) 


V 


i o 


3 This is in consonance with the mass term being equivalent to a helicity flip. 
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with o i being the first Pauli matrix. We will approximate Eq. m by a differential equation 
by first writing 


( 4>+(z-At,t-At)\ 

= Ttj}{z,t), 

+ At ,t — At) J 

with the “transfer” operator T given by 


( 5 ) 



( 6 ) 


where <73 = diag(l, —1) is the third Pauli matrix. Then the difference equation Eq. ((21) 
takes the form 

= TTil)(z,t). (7) 


The difference equation becomes a differential equation if we limit ourselves to the zeroth 
order (given by the identity ij)(z, t ) = ip(z, t)) and the first order term in At. The first order 
equation is 


d d 

z— t) = m<Tii(}(z, t ) - %a z —^{z, t ), 


( 8 ) 


the Dirac equation 4 in (1 + 1) dimensions! There is no spin in the SG^l, 1) case, the little 
group of being trivial. I 11 passing we mention that, for the Dirac equation to hold to 
all orders in At, due to the fact that oq and 03 do not commute with each other, it is 


necessary to replace tFT ii 




will not pursue this issue any further and will assume that At is sufficiently small that higher 
order terms are negligible, i.e., the Dirac equation is a good approximation to the original 
difference equation. 


III. HIGHER SPATIAL DIMENSIONS 

We start by reminding ourselves that, for the Dirac equation, the velocity operators 

-k[x, H] = ca not only have eigenvalues ±c, but they also anticommute with each other . 5 

The latter fact makes the generalization of our approach to more than one spatial dimension 

4 In our representation of the Dirac matrices, the positive-energy spinor for the plane-wave solution takes 
the form, aside from a normalization constant, u(p ) ~ (—ay \fp 2 + m 2 — i<J 2 P — m)u(0), with u( 0) being 
the two-component spinor (1, —1), and the negative-energy spinor v(p) ~ (tri \Jp 2 + to 2 + ioqp — m)v{ 0), 
with u(0) = (1,1). 

5 Here an analogy with spin can be made. Just as one cannot specify all three components of the spin 
simultaneously in quantum mechanics without running into inconsistencies (of predicting a spin V3 times 
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non-trivial. Before we proceed to the (3 + l)-dimensional case, let us first discuss the two 
spatial dimensional SO (2,1) case. 

As in the preceeding section (for the SO(l, 1) case), we use a 3 to give the dependence 
on the z coordinate and 0 \ for the flip operator. Of the three Pauli matrices, we have only 
cr 2 left; so let us call the second spatial coordinate the y coordinate. Now the problem is to 
express the dependence of ip on y in the sense that a 3 gives the dependence on z. 

To solve this problem we appeal to rotational invariance and make explicit use of the 
spin-| property of ip. Let 1Z be the rotation over 7 t/2 from the y axis to the z axis. Then 
U{71), which represents this rotation, is given by U(1Z) = 4=(1 iai). Thus, to the cr 3 Jl 
term in the “transfer” operator T we add 


d 


U~\n)a 3 —U(n) = a 2 


d_ 


and obtain 


.d ( d d \ . . 

l ~dr^ V ' ^ = \ ai ~ m2 ~dy ~ W3 'd z I yj y ’ Z ’ f ' 


(9) 


( 10 ) 


Note that the flip operator, which is used in the preceeding section to invert the ^-motion, 


also inverts the y-motion. In Dirac’s notation[lj, we identify (7\ = a mi cr 2 = cx 2 , a 3 = a 3 . 

Now we recall the general rulejsj that spinors in 2 n dimensions and in 2n+ 1 dimensions 
have 2 n components. Thus for the case of three spatial dimensions and one temporal dimen¬ 
sion (50(3,1)), we need 4-spinors. And we need, besides a m , a 2 , a 3l an extra a, all four 


of which anticommute with one another and have eigenvalues ±1. Following Dirac [jj, we 
introduce two independent sets of Pauli matrices p and a. The p’s and the cr’s anticommute 
among each set, whereas the p’s and the cr’s commute. As we want to make a m look like a 
flip matrix, we pick 

^OOl' 


Pl&l 


(ii) 


0010 
0100 
1000 

We complete the set (by the same argument we have used above for the 50(2,1) case) with 


V 


/ 


01.1 — P 2 V 1 , 0(2 — p 3 <Ji, oi 3 — cr 3 l. 


( 12 ) 


bigger in some directions), one cannot specify all three components of the instantaneous velocity simul¬ 
taneously without running into inconsistencies of predicting a tachyonic speed of \/3 times the speed of 
light. 
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Here 1 is the 2x2 unit matrix. In passing, we mention that it is easy to treat the case of 
(4 + l)-dimensional spacetime, as we can now identify ol\ with o^l. 

For the case of (3 + 1) dimensions, the equation we obtain is 





d d d\ , 

-- at 3 i— ip. 

ox oy oz) 


(13) 


It is of relevance to remark that the last three terms on the right hand side of Eq. m 
approximate the small finite steps of motion with the speed of light before the event (x, y, z, t ) 
is reached. The term a m m represents the unitary transformation e~ lFmdt which takes place 
at that event (see Eqs. © and ©)• Thus the electron is not aware of the fact that it has a 
mass until it hits a lattice site. If it has no mass, then it is not flipped and it moves at the 
constant speed of light. If it is massive and is at “rest”, then it must be that the electron 
zigzags around with the speed of light and returns to its original spatial lattice site and 
wanders around again and returns again etc. 


IV. INTERACTIONS WITH EXTERNAL FIELDS 

To put the Dirac equation in a covariant form, we follow the usual procedure of writing 
a m = 7 ° (with ( 7 0 ) 2 = 1) and multiplying Eq. (1T31) by 7 0 to yield 

m)ip = 0, (14) 

where (/i running over 0,1, 2, 3) 

7 m = (15) 

with a a = /, the identity matrix and 0„ = (|. £ £, &). FVom the way we have derived 
the Dirac equation, we can trace the term to the “transfer” of the electron at the 

speed of light between the lattice sites while the m term comes from the “flip” unitary 
transformation at the lattice sites. 

The introduction of an electromagnetic held A^ is straightforward by using the prescrip¬ 
tion of minimal substitution in Eq. (ED 

id iL =7- id fl + eA^. (16) 

Although the ea^A^ term goes together with the term in the minimal subsitution rule, 
it is tempting to keep the d^ term identified with the transfer between lattice sites and put 
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the ea^A^ term together with the Fm flip term as taking place at the lattice sites. (But 
we should keep in mind that, since the Dirac matrices do not commute among themselves, 
beyond the first order term, there is a difference between associating the interaction term 
with the “transfer” operator T or the “flip” operator JF. 6 ) 

To incorporate gravitational interactions one needs the tetrad (or vierbein) formalism^. 
One introduces at every event x a set of local inertial coordinates with a tetrad e' J : of 

Q 

axes labelled by the Minkowski index a, b, c running over 0,1, 2, 3. 7J Then the metric in 
any general noninertial coordinate system is given by g tw = e a ^g a be b u in terms of the flat 
Minkowski metric 7] ab . Gravitational interactions are introduced via the substitution rule 

rd a =>re u a (d u -Ucu bcu a bc ), (17) 

where cu b u = (d u e£ + T() A e^)e^ in terms of the affine connection T^ x and a bc = \i [y 6 , y c ]. At 
every spacetime lattice site labelled by x, we have a tetrad e%. In our interpretation, the 
electron travels with the speed of light between lattice sites; this is represented by 
Then at the lattice site there is a unitary transformation which, in addition to the mass 
“flip”, now contains the interaction term 7 “ e^iVbcuO' bc ■ 

V. DISCUSSIONS 

We have presented a novel derivation of the Dirac equation, hoping to shed new light 
on the physics of the electron. Motivated by the distinct possibility that the underlying 
spacetime is discrete at small scales, we have started with a discrete “cubic” lattice. The 
resulting Dirac equation emerges as the lowest nontrivial order of approximation. Thus the 
observed Lorentz invariance does not preclude the existence of a discrete spacetime at small 
scales. 

Is our approach useful? We think so. (1) The very fact that the underlying spacetime is 
discrete means that there is automatically an ultraviolet cutoff which may be used to ame¬ 
liorate divergence problems in nonrenormalizable theories like (perhaps) quantum gravity. 
(2) Our starting point is a difference equation rather than a differential equation. While dif¬ 
ference equations are more tedious to deal with analytically, they may hold some advantages 
in numerical calculations. 

6 For the 50(1,1) case, the A M term can be incorporated into the flip operator in Eq. ©. 


7 



We conclude with some speculations and a couple of open questions. In the sce¬ 
nario we have proposed, the electron travels between lattice sites with the speed of 
light. Only at the lattice sites does the electron “feel” its mass and perhaps also the 
presence of all external fields. 7 (Since it is a Yukawa-type interaction which, via the 
Higgs mechanism, generates mass for the electron, it seems reasonable to assume that 
at least Yukawa-type interactions take place only at the lattice sites where the mass 
operator makes its presence felt.) But if gravitational interactions also take place mainly 
at the lattice sites, does that mean spacetime vertices somehow play an important role 
in concentrating curvature? And if so, how is this description of geometry and topology 
related to the Regge calculus jsj], for example? These problems deserve further investigations. 
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